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Abstract 



Series expansions at T = are used to study properties of the half-filled 
t-J ladder doped with one or two holes, and at quarter filling. Dispersion 
curves are obtained for one-hole symmetric and antisymmetric (bonding and 
antibonding) excitations and for the two-hole bound state. The line in the 
phase diagram that separates bound and unbound states is determined. For 
quarter filling we compute the ground state energy and estimate the location 
of the phase separation line. Comparisons with other numerical and analytical 
results are presented. 

PACS Indices: 71.27.+a, 71.10Fd 
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I. INTRODUCTION 



The discovery of materials containing S = ^ ions which form a 2-leg ladder structure 
has stimulated a great deal of theoretical work to understand the unique and surprising 
features of these systems (for reviews see ) • Early work on magnetic properties confirmed 
the existence of a spin-gap in the 2-leg ladder with antiferromagnetic (AF) interactions, and 
hence a qualitative difference between even and odd legged spin ladders. More exotic spin 
ladders, with frustration, dimerization and other features have been studied more recently. 



Even more interesting behaviour may be expected if the system includes charge degrees 
of freedom. This can be achieved by doping, to create a system of strongly correlated 
mobile holes, as in the cuprate superconductors. Indeed an understanding of this (arguably 
simpler) system may yield important consequences for a better understanding of high 
superconductivity. The simplest model, in this context, is the t-J model away from half- 
filling. Considerable work has been done on the t-J model on a 2-leg ladder, through exact 
diagonalizations |[T0|-[T6|, the density matrix renormalization group approach and 
approximate analytic theories |T^-^ . A number of broad features have emerged from this 
work. 

The t-J ladder is shown in Figure 1. The Hamiltonian is 

H = J ^(Sj^a ■ Sj+i,a — -ni^ani+l^a) + J 1. '^{^i,! " ^i,2 — j'^i, 1^^1,2) 
i,a i 

-t E P{cla,.Ci+i,a,a + H.C.)P - t^Y. ^(4l,<xC,,2,. + //.C.)P (l) 
i,a,a 1,(7 

where i runs over rungs, cr (=| or |) and a (=1,2) are spin and leg indices, c|^^ (^,0,0-) is the 
creation (annihilation) operator of an electron with spin a at the ith site of the ath chain, 
and Si^a = \(^i^a,a^a,pCi,a,li deuotes the 5* = I spin operator at the ith site of the ath chain. 
P is a projection operator that ensures that doubly occupied states are excluded. 

To understand the nature of the ground state and excitations of the system it is helpful 
to consider the limiting case J±,t± ^ J,t when the states on the rungs of the ladder form an 
appropriate basis. Many of these features persist even to the isotropic case J± = J, t± = t. 
At half-filling the charge degrees of freedom are frozen out and the system is equivalent to 
the Heisenberg spin ladder. The ground state is a gapped spin liquid with each rung in 
a spin singlet state. Spin excitations originate from exciting one or more of the rungs to 
triplet states. 

When the system is doped, leading to removal of electrons, all nine states of a rung must 



be included. These are discussed in Ref. [11|, and for completeness are illustrated in Table I. 



A single hole (of either spin) will go into a bonding or antibonding state on a rung and this 



state can propagate along the ladder. Troyer et al. [|Tl|] have obtained one-hole excitation 
spectra from exact diagonalizations of systems up to 20 sites. They argue that the free spin 
created remains bound to the hole and the resultant quasiparticle carries both spin and 
charge. This is unlike the situation in the t-J chain where spin-charge separation occurs, in 
accordance with Luttinger liquid theory. 

Many authors have considered two-hole states. It is energetically favourable for the two 
holes to form a bound state, on the same rung. The lowest two-hole excitation arises from 
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a coherent propagation of the two holes along the ladder. This excitation carries charge 2e 
and no spin. There are other, more complex, excitations which are described by Troyer et 
al. Q. 

The nature of the system at finite doping is also of considerable interest. Poilblanc et 



al. [0 have conjectured a T = phase diagram in the plane of J/t vs electron density 
n. For large J one expects phase separation, where the system separates into hole-free and 
hole-rich regions. For smaller J/t and at or near half-filling, there is a "CISO" phase p2 



with one gapless charge mode and no gapless spin modes. The spin-gap will decrease, as n 
is decreased from 1, and a transition to a Luttinger liquid phase is predicted. At small n 
electron-pairing rather than hole-pairing will occur. Miiller and Rice have extended this 
phase diagram through the inclusion of a Nagaoka phase at low doping and small J, and 
a "C2S2" phase with gapless charge and spin excitations at larger J, between the Nagaoka 
and spin-gapped phases. Much of this remains somewhat speculative. 

In this paper we study the t-J ladder via the method of series expansions. This approach 
is complementary to other numerical methods and is able to provide ground state properties, 
quantum critical points, and excitation spectra to high accuracy. Its other advantage is that 
one deals, at once, with a large system and hence there are no finite size corrections needed. 



We have used this approach recently in a study of the t-J model on the square lattice, 
and we refer to that paper for details of the method and references to previous works. Early 
perturbative studies of the t-J model only considered the lowest order terms. The 
emphasis of our work is on one- and two-hole states, and we display excitation spectra for 
these. The method is not well suited to handle finite doping and hence we are not able to 
add to knowledge of the phase diagram for general n. We do however consider the special 
case of quarter-filling, and obtain an estimate of the boundary of phase-separation in this 
case. Wherever possible we compare our results with previous work. 

This paper is organized as follows. In Sec. II, we consider the system with one and two 
holes. In Sec. Ill, we study the system at quarter-filling. The last section is devoted to 
discussion and conclusions. 



II. ONE-HOLE AND TWO-HOLE STATES 

The series expansion method is based on a linked cluster formulation of standard 
Rayleigh-Schrodinger perturbation theory, with the Hamiltonian (|l]) written in the usual 
manner as 

H = Ho + xV (2) 

where Ho is a solved, unperturbed system and xV represents the remaining terms. The 
most obvious choice is to include all of the rung terms in Ho 

Ho = J± ■ S':2 - ^ni^i^i^a) ~ ^-L H ^(4l,aCi,2,<7 + H.c.)P (3) 

i i,a 

and the coupling between rungs as V 

xV = J ■ Si+l,a - ^ni^ani+l,a) " t ^ P(c^„Ci+i,a,a + H.C.)P (4) 

i.a i.a.a 
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As discussed above, the eigenstates of Hq are direct products constructed from the nine 
possible rung states, shown in Table I. To compute the perturbation series we fix the values 
of t±/ J± and t/ J and derive an expansion in powers of x = J/ J± for the quantity of interest. 
The series is then evaluated (summed) at the desired value of J/ J± by using standard Fade 
approximants and integrated differential approximants We set J_l = 1 to define the 

energy scale. To improve the convergence of the expansion we also try the effect of adding 
a term to Ho and subtracting it from V, for example 

K = {J± + r) 5](S,- 1 ■ S,,2 - ^n,,in,,2) - t± E Piil,.c^,2,a + H.c.)P (5) 

i i,a 

xV' = xV - r E(Si,i ■ Si,2 - ^^i,i?^i,2) (6) 

i 

and adjusting r to obtain best convergence. Other types of term could be added and 
subtracted in this way. 

A slightly different formulation takes only the rung spin Hamiltonian as the unperturbed 
part 

^^0 = J± Yli^iA ■ - ^^i,l^i,2) (7) 

i 

with all of the remaining terms, including the rung hopping, included in the perturbation. 
Both methods have been used, as appropriate. 

The series method can be used to compute the ground state energy and other ground 
state properties such as magnetization, susceptibility, and correlations. In the present study 
we are mainly interested in the nature of excitations. Rather than calculate the ground 



state energy we have used the method of Gelfand to compute the excitation energy 
directly. For the one-hole excitations we start with an initial state with the hole in either 
the bonding or antibonding state on a particular rung and all other rungs in spin singlet 
states. The perturbation will then mix this with other states in which the hole moves along 
the ladder within a fluctuating spin background. For the two-hole case the initial state has 
both holes on the same rung. During the evolution of the system these holes may separate. 
Those final states which have the holes together on the same final rung contribute to the 
two-hole excitation spectrum. The series have been computed to order (J/J±)^^ for the one 
hole case and to order {J/J±y^ for the two-hole case. In special cases we have obtained 
considerably longer series. In the case of static holes (t = 0) we have extended the series to 
order (J/J±)^^ for both one- and two-hole cases, while for J = we have obtained one-hole 
series to (t/J±y^ and two-hole series to (t/J±y^. 

The series are available on request. For the interest of the reader we give some of the 
low order terms for both one and two-hole series in Appendix A. As usual with series meth- 
ods the convergence is judged by consistency of different Fade (or integrated differential) 
approximants. In the present case convergence worsens for increasing J/J± and t/J±. Error 
estimates are shown in the figures displaying our results. 

Figure 2 shows the excitation spectra for one-hole symmetric and antisymmetric states for 
J/ J± = 0, 0.5, 1.0 and various t/ J±. To avoid excessive data we have set t = t± throughout. 
The quantities graphed are the gaps AE^^ / J±, AE^^/J±, where for instance 
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AE^^ = E^^ - Eo (8) 

and Eq is the energy of the half-filled ground state. We draw attention to the following 
features. 

• The series yield smooth dispersion curves throughout the 1-d Brillouin zone. There 
is consistency between different approximants and we believe the error is small. The 
only exception is for some of the antisymmetric states where the figure shows larger 
error bars. 

• For fixed J/ J± the energies decrease with increasing t. The kinetic energy gain, pro- 
portional to t, exceeds the loss in potential energy from breaking antiferromagnetic 
bonds. The bandwidth also increases with increasing t, as expected. The gap vanishes 
or goes negative when t/Jj_ gets large, i.e. the 1-hole state lies lower in energy than 
the half-filled ground state. 

• The main effect of increasing J is in the position of the minimum. For the symmetric 
states this is at A; = tt for small J, and moves to a value near = 27r/3 for J greater 
than a critical value around J/J± — 0.6, almost independent of t. 

• The antisymmetric state shows greater sensitivity to the value of J, with a qualitative 
change in the one- hole dispersion curve between J = and 0.5. 



These dispersion curves are in broad qualitative agreement with those obtained by other 
methods [pT] , ^ , but with significant quantitative differences. In Figure 3 we reproduce 
our results for the one-hole symmetric state spectrum for J/J± = 1, and compare with the 



diagonalizations of Troyer et al. |TT|, and an approximate analytic result of Sushkov pT 
The dispersion curve obtained from exact diagonalization has a very similar shape to our 
one, but lies at somewhat lower energies. The difference may be due to finite size effects in 
the diagonalization results. The analytic calculation, on the other hand, lies at high energies 
and appears to overestimate the dispersion. The discrepancy is greatest at = 0. Note that 



ref. I^T)] does not include the —nirij/A terms in the t-J Hamiltonian and hence a constant 
0.75 must be added. 

We now turn to the two-hole excitations. Figure 4 shows the excitation spectra, again 
for J/J± = 0,0.5,1.0 and for various values of t/J±. The quantity graphed is again the 
2-hole gap 

where E'^^ is the energy of the lowest-lying singlet 2-hole state. The following features are 
apparent 



• The minimum now occurs at A; = for all t, J. 

• The curves are quite flat for small t/J±, as expected, and develop greater dispersion 
and more structure as t increases. The matrix element for two holes to hop from one 
rung to the next is of order t^. 
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• The energy at k = ir depends only weakly on t. For the special case J = the energy 
at = TT is in fact completely independent of t. It is possible to prove this results 
quite simply by looking at matrix elements of two-hole hopping processes. We do not 
have a simple physical explanation. 

• The two-hole gap vanishes or goes negative for small kx when t/J± gets large. 
The two-hole binding energy can be computed from our excitation energies by 

Eb = AE^'' - 2AE^'' = E^'' - 2E^^ + E^ (10) 

where AE^^ is the lowest 1-hole energy, which is the symmetric one in our case. Negative 
Ef, corresponds to binding of the two holes. Figure 5 shows Eh versus t/ Jj, for various values 
of J/ J±. For t = all the curves converge to the vicinity of —1, which is the exact (trivial) 
binding energy for J = 0. For larger J this value increases slightly (to approximately -0.87 
for large J), but this is below the resolution of this figure. For small J the binding energy 
decreases in magnitude for increasing t, and the two holes no longer bind beyond a critical 
tc, which is itself an increasing function of J. For larger J binding seems to occur for any 



t. Our results are numerically consistent with those of Hayward et al. |14] and Troyer et al. 



11 . For the isotropic case J = J±, the two hole state is bound for all values of t. 



III. THE QUARTER-FILLED CASE 

We have obtained a number of results for the case of quarter filling, i.e., where there 
is one electron (or hole) per rung. This value of electron density lies within the Luttinger 
liquid (ClSl) phase in the conjectured phase diagram ]12|], but rather near and parallel to 
the line separating this phase from the spin-gapped phase. There has been recent interest 



in this case as a model for NaV205 27-2S 



To develop the perturbation expansions for quarter-filling we start from an unperturbed 
state which is a direct product of alternating spin-up and spin-down rung bonding states. In 
the absence of V this is but one of 2^/^ degenerate states, because the spin-up and spin-down 
rung bonding states have the same energy. To remove this degeneracy, we have considered 
the effect of the perturbation \^ on a system of two rungs, and have included the diagonal 
terms in the unperturbed Hamiltonian. This splits the degeneracy and favours a Neel type 
ordering along the ladder. To improve the convergence of the expansion we also try the 
effect of adding the following term AH to Hq and subtracting it from the perturbation V, 

AH = r'j:i~iysum - \sLm (n) 

i 

where i runs over rungs, and S^^^^{i) is the total spin in 2;-direction for rung i. We adjust the 
coefficient r' to obtain best convergence. The series have been computed to order x^^ for 
the ground state energy. 

In Figure 6 we present results for the ground state energy as a function of J/ J± (taking, 
as usual, t± = t). Also shown is the curve for the half-filled ground state energy (divided by 
a factor of 2). The crossing points provides an estimate of the locus of the phase separation 
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line, where it becomes energetically favourable for the electrons to cluster together. These 
results are plotted and discussed in the following section. 

An interesting limiting case is where ^ t, J, J±. In this limit there is a mapping to 
a single = | antiferromagnetic Heisenberg chain with Jgfj = J/2. This mapping arises 
from replacing the rungs with effective S = \ spins, and can be established by a detailed 
investigation of matrix elements. The corresponding effective Hamiltonian for this mapping 
is 

H'^ = -t^N/2 + Jeff E(S« ■ S.+i - 1/4) (12) 

i 

In Figure 7 we show the ground-state energy per spin versus J/t± for two choices of 
{t, J±), together with the energy per site of the corresponding chain 

Eq/N = -t±/2 - J In 2/4 (13) 

The agreement is exact near J = and quite good even up to J/t± = 1 provided t and J± 
are small. 



IV. DISCUSSION AND CONCLUSIONS 

We have carried out the first comprehensive series study of the t-J ladder system, which 
is applicable to real materials as well as being of theoretical interest in its own right. We 
have considered the system near half-filling (one and two holes) and for exactly quarter 
filling. Our results provide a comparison for other numerical and analytical studies of the 
t-J ladder, especially exact diagonalizations for small systems which have been extensively 
used. 

Both one-hole and two-hole states propagate as quasi-particles with well defined excita- 
tion energies, and our estimates of these are, we believe, the most accurate available. They 
are in qualitative agreement with exact diagonalization results but the latter are for small 
systems and may have substantial finite size correction terms. We have obtained precise 
numerical estimates for the two-hole binding energy. 

a) Half-filling 

Figure 8 shows the resulting "phase diagram" for the half-filled case. From Figure 5 we 
can estimate the critical line at which the 2-hole binding energy vanishes, shown as line (A) 
in Figure 8. Below the line two holes bind, above the line they are unbound. Also shown 
are the lines where the 1-hole gap vanishes (B), and the 2-hole gap vanishes (C). 

The binding of two holes is a necessary requirement for phase separation, where holes 
clump together to form separate hole-rich and hole-poor regions; but it is not sufficient. The 



usual criterion for the phase separation boundary is taken to be |1TT|JT^ the point where the 
inverse compressibility vanishes: 

^ ^ [E{N + 2) - 2E{N) + E{N - 2)] (14) 



where N is the number of electrons (or holes). In other words, phase separation occurs when 
the energy surface becomes locally convex with respect to the addition or subtraction of 
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hole pairs; whereas the binding energy ([TOD measures the convexity with respect to addition 
of single holes. We are unable to compute the quantity ([T^) . Tsunetsugu et al. [|ll|] and 
Hayward et al. use criterion ([1^ to find the phase separation boundary at line (D) in 



Figure 8, very much lower than the 2- hole binding boundary (A). 



This model is also believed |12,11] to exhibit a "CISO" phase above the phase separation 



boundary with gapped spin excitations but gapless charge excitations. Troyer et al. |Tl 
discuss how the bound hole pairs form a gapless band of charge fiuctuations, and how the 
dispersion relation for the pair becomes linear near k = 0. The most naive criterion for the 
boundary of the CISO phase would be the vanishing of the two-hole gap, line (C) in Figure 
8. That would imply the possibility of a narrow intermediate region between lines (C) and 
(D), corresponding to a "COSO" phase where phase separation has not yet occurred. Other 
authors define the "charge gap" 

A, = \[E,{N^ + 1, N{) - 2EoiN^, iVj + E^iN^ - 1, iVj] (15) 

which cannot be computed with our current techniques of series expansions, 
b) Quarter-filling 

From Figure 6, we can estimate the boundary at which 

Eo{n = i) = Eoin = l)/2 (16) 

where n is the number of electrons per site (n = 1 for half-filling, n = 1/2 for quarter- 
filling). This is a sufficient but not necessary criterion for phase separation, implying a 
'global" convexity of the energy surface, so that the system could save energy by separating 
completely into a hole-empty and hole-full region. Estimates based on this criterion are 
given by the dotted line in Fig. 9. Estimates based on the inverse compressibility |TT 



are 



shown by the open points in Figure 9: they lie a little above our estimates, but quite close 
to it, and are hardly distinguishable within present accuracy limits. 

We have noted on a mapping (exact at the limit t, J, J± <^ t±) from the quarter-filled t-J 
ladder to an S* = | antiferromagnetic Heisenberg spin chain. This idea can be generalized to 
other systems, such as the dimerized t-J chain or dimerized t-J model on a square lattice, 
and we will report on this in future work. 

Note added in proofs: After submission we learnt of recent work of Bose and Gayen 



in which they obtain exact results for 1 and 2-hole states in a t — J ladder with additional 
diagonal interactions. The physics of two systems appears similar. 
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APPENDIX A: 



If we take the term in Eq. (|^) as the unperturbed Hamiltonian, and all the remaining 
terms as the perturbation, we can get a double series for the excitation energies. The series 
up to order fj^ {i+ j < 4) for the excitation spectra of the one-hole symmetric state AE^^ 
is (with J_L = 1 and t± = t) 

AE^^ = 1 + (J/2 -t) + (3JV8 - 3tV2) + (3JV16 + 3tJ^/A - 9tV/4 + 3t^) 
+ (-3J732 + 3tJ^/A - 69tVV32 + 9t^J - 105tVl6) 

+ {t - 3t J/2 - 9t JVS + 3tV - 9tV4 + 9t JV32 + 9t'^J^/2 - 15tV/2 + 9t^) cos{h) 
+ (15tVVl6 + 3tV/2 - 3t^/2) cos(2A;^) + ■ ■ • (Al) 

The series for the one-hole antisymmetric state is related to that for the symmetric case 

by 

AE^^{t, J, k^) = AE^^{-t, J, TT - k^) (A2) 

The two-hole excitation energy AE"^^ is 

AE"^^ = 1 + J + 3 JV4 - 4t2 + 3 JVS - 3 JV32 + 3t^J^/8 + 

+ (-4^2 + 3t2 jV4 + 6t^) cos(A;^) + llt^ cos(2A;^)/2 + • • • (A3) 
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FIGURES 



FIG. 1. Stucture and parameters of the t-J ladder. 



FIG. 2. Excitation spectra for one-hole symmetric and antisymmetric states for 
J/J± = 0,0.5,1.0 and various t/J±. In each case the dispersionless energy for t = is also 
shown. 



FIG. 3. Comparison of our results (points connected by dashed line) for the one- hole symmetric 
spectrum with exact diagonalizations (Troyer et al. |11|, solid curve) and an approximate analytic 
theory (Sushkov pl|, dotted curves). 



FIG. 4. Excitation spectra for two-hole states for J/ J± = 0,0.5, 1.0 and various t/J±. In each 
case the dispersionless energy for t = is also shown. 



FIG. 5. Binding energy of two holes vs t/J±. The J = case shows estimates for different 
integrated differential approximants, which begin to spread out for t/J± ^ 1. The points where 
Eh = give the two-hole binding line shown in Fig. ^. Also shown are the results of Hayward et 
al. 1 14] for the isotropic case J/J± = l (the big full points connected by the dot-long dash line). 



FIG. 6. Ground state energy per site for quarter-filling versus J/J± for various t/J±. The solid 
line is the ground state energy for the phase separated state. The crossing points yield the phase 
separation line shown in Fig. ^. 



FIG. 7. Comparison of the ground state energy per site for quarter-filling for the case of t = 
(note that the results are independent of J± for a small value of J±) and the case of J± = t = J 
with the results of an approximate mapping to a Heisenberg spin chain (solid line). 



FIG. 8. Phases and critical lines for the t-J ladder at half-filling, in the plane of t/J± versus 
J/J± (see text). Two holes do not bind in the region marked NB. 



FIG. 9. Phase separation boundary at quarter- filling in the t/J± versus J/ J± plane. The region 
marked PS is where phase separation occurs. The open points mark the boundary value found by 
Hayward and Poiblanc ||ll[ for periodic (PBC) and open (OPC) boundary conditions on 8-rung 
finite chains. 
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TABLES 



TABLE I. The nine rung states and their energies. 



No. 


Eigenstate 


Eigenvalue 


Name 


1 


;fe(in>- lit)) 


-J± 


singlet 


2 


lii) 





triplet = -1) 


o 
.) 


1 / ht 1 \ 1 1 1 t \ ^ 

lif)) 


U 


tdipict. (..S|„| = {)) 


4 


ITT) 





triplet = 1) 


5 


1 00) 





hole-pair singlet 


6 


7l(|Oi)+|iO)) 


-t± 


electron-hole bonding (S^ot — 


7 


7i(|0T)+IT0)) 


-t± 


electron- hole bonding (S'fot = ^) 


8 


;^(|Oi)-|iO)) 


t± 


electron-hole antibonding {S^^^ ~ ~l) 


9 


^(loT)-ITO)) 


t± 


electron-hole antibonding (^tot = |) 
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